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Abstract 


In this study, we present a catalog of molecular clumps extracted from '*CO (J = 1 — 0) emission data of the Milky 
Way Imaging Scroll Painting (MWISP) project. The data covers the inner Milky Way within the longitude range 
10° <1 < 20? and the latitude strip of |b| < 5°25. The workflow for the extraction of clumps, namely Facet-SS-3D- 
Clump, consists of two parts: the identification of clump candidates and their verification. First, Facet-SS-3D- 
Clump employs FacetClumps to identify clump candidates. Subsequently, high-confidence clumps are obtained by 
cross-matching with the clumps detected by other algorithms, such as dendrogram. Second, these high-confidence 
clumps are used as prior knowledge to train a semi-supervised deep clustering approach, SS-3D-Clump, which is 
applied to verify clump candidates detected by FacetClumps, providing confidence levels for the molecular 
clumps. Finally, the catalog comprising 18,757 molecular clumps was obtained using Facet-SS-3D-Clump, and the 
catalog is 90% complete above 37 K km s™!. We observe a significant deviation of the mean Galactic latitude for 
clumps within |b| < 2° from the midplane, with b = — 0°110. We found that 82.3% of the dust clumps correspond 
to PCO clumps by matching with Herschel infrared dust clumps. In the future, Facet-SS-3D-Clump will be applied 
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to detect ^CO clumps in the entire MWISP data. 
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1. Introduction 


It has been accepted that stars are formed in molecular 
clouds (MCs; e.g., Shu et al. 1987; Krumholz & McKee 2005; 
Zinnecker & Yorke 2007; Krumholz et al. 2009; Nejad- 
Asghar 2011; Kónyves et al. 2015; Basu 2016), while the 
molecular clumps that condensed from the ambient cloud 
represent a critical step in the star formation process (Hacar 
et al. 2013). The molecular clump formation also seems to play 
an essential role in determining the mass of the final star since 
the distribution of masses among starless cores in a cloud 
mimics the stellar initial mass function (IMF; e.g., Motte et al. 
1998; Alves et al. 2007; Nutter & Ward-Thompson 2007; 
Kónyves et al. 2010; Benedettini et al. 2018; Bresnahan et al. 
2018; Arzoumanian et al. 2019). Although great progress has 
been made in determining the form of the IMF (e.g. 
Chabrier 2003; Nutter & Ward-Thompson 2007; Olmi et al. 
2009; Kroupa et al. 2013; Kónyves et al. 2015), a detailed 
explanation of the form and possible environment depends on 
understanding the nature and evolution of the molecular clumps 
(Marsh et al. 2016), especially dense molecular clumps. Thus, 
the clumps of molecular interstellar medium (ISM) surveys 
play an important role in understanding the location and mode 
of star formation in the Milky Way. Since the discovery of CO 
molecules in the 1970s, numerous survey projects have been 


carried out successively, such as the Galactic Ring Survey 
(GRS; Jackson et al. 2006), the Exeter-FCRAO CO Galactic 
Plane Survey (Mottram & Brunt 2010), the Milky Way 
Imaging Scroll Painting (MWISP; Jiang & Li 2013), the Three- 
mm Ultimate Mopra Milky Way Survey (ThrUMMS; Barnes 
et al. 2015), the CO Heterodyne Inner Milky Way Plane Survey 
(CHIMPS; Rigby et al. 2016), the Structure, Excitation, and 
Dynamics of the Inner Galactic Inter-Stellar Medium (SEDI- 
GISM; Schuller et al. 2017) survey, the Forgotten Quadrant 
Survey (FQS; Benedettini et al. 2017), etc. 

All of those surveys, with different sensitivity and spatial 
resolution, could help us to detect the distribution of the 
molecular component with different scales (Benedettini et al. 
2021): from dense clouds and filamentary structures (André 
et al. 2014) to pre-stellar clumps and young stellar objects 
(YSOs). The MWISP (e.g., Zhang et al. 2014; Sun et al. 2015; 
Du et al. 2017; Sun et al. 2017; Su et al. 2019) project 
implemented with the Purple Mountain Observatory Delingha 
(PMODLH) 13.7 m telescope (Zuo et al. 2011) is dedicated to a 
new large-scale survey of molecular gas in PCO (J=1-— 0), 
BCO (J—1— 0), and C^O (J=1-— 0), which targets the 
northern Galactic Plane within —10°25 «/« 250725 and 
—5°25 < b < 5°25 and several other regions of interest. The 
goal is to develop our knowledge of molecular gas's precise 
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content and distribution in the Milky Way, the physical and 
chemical processes converting tenuous interstellar gas to dense 
molecular gas, and the rules governing the formation of stars 
from MCs. 

With the progress of sky survey projects, unsupervised 
machine learning-based algorithms for detecting molecular 
clumps have also emerged (e.g., Rosolowsky et al. 2008; 
Berry 2015; Luo et al. 2022; Jiang et al. 2023). The 
dendrogram algorithm, developed by Rosolowsky et al. 
(2008), is well-suited for illustrating changes in the hierarchical 
structure of isosurfaces within molecular line data cubes as 
contour levels vary (Rani et al. 2023). It has been widely used 
in continuum, atomic hydrogen, and molecular line data 
(Cheng et al. 2018; Takekoshi et al. 2019; Nakanishi et al. 
2020; Zhang et al. 2021). Luo et al. (2022) utilized the local 
density clustering method, inspired by the concept from Alex 
Rodriguez (2014), to identify clumps as local dense regions 
embedded within molecular gas with lower average bulk 
density, as suggested in (Blitz & Stark 1986; Lada 1992; 
Bergin & Tafalla 2007). Additionally, a Multiple Gaussian 
Model (MGM) is utilized to simultaneously fit overlapping 
clumps, allowing for the extraction of clump parameters. Jiang 
et al. (2023) presented FacetClumps, a Gaussian facet-based 
model designed for fitting local surfaces. The extremum 
determination theorem of multivariate functions is applied to 
determine clump centers. Then, local regions of the centers are 
clustered to identify clumps by considering connectivity and 
minimum distance. The experiments (Jiang et al. 2023) show 
that FacetClumps exhibits excellent recall and precision rates. 

Accumulating survey data and developing automatic detec- 
tion algorithms for molecular clumps have made it possible to 
conduct a comprehensive census of these clumps. For example, 
Wu et al. (2012) expanded the horizon of cold astronomy by 
studying the Planck cold clumps. Marsh et al. (2016) published 
a catalog of dense cores of the Taurus star-forming region 
derived from Herschel SPIRE and PACS observations. 
Benedettini et al. (2020) presented a catalog of MCs extracted 
from the FQS CO (J = 1 — 0) spectral cubes. Typically, 
researchers will mitigate false positives through visual inspec- 
tion, such as Rigby et al. (2019) visually inspecting each 
CHIMPS clump by three independent reviewers and assigning 
a reliability flag. Rojas et al. (2022) visually inspected and 
classified the lens candidates into two catalogs using high- 
resolution imaging and spectroscopy. At the same time, 
researchers need to devise meticulous strategies to eliminate 
potential subjective factors introduced during the visual 
inspection process, which will require a significant amount of 
work and time. Therefore, an automated verification method as 
a substitute for manual verification becomes increasingly 
necessary in the era of extensive data. Luo et al. (2024) 
propose a semi-supervised deep clustering method for 
molecular clump verification, namely SS-3D-Clump, which 
extracts deep futures of clumps and classifies these features 
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through a clustering algorithm to obtain a pseudo label. 
Subsequently, it uses these pseudo labels as supervision to 
update the weights of the entire network. It can leverage 
unlabeled samples through semi-supervised learning to 
enhance the generalization ability of SS-3D-Clump and has 
shown remarkable performance. 

In this paper, we present a comprehensive workflow for 
detecting and verifying OD clumps in the MWISP project. 
The workflow, namely Facet-SS-3D-Clump, first employs 
FacetClumps to obtain candidates for molecular clumps and 
then uses SS-3D-Clump for their verification. To validate the 
reliability of the workflow, we selected a subregion within 
10°<1<20° & —5°25<b<5°25 from the MWISP as 
experimental data, which includes some active high-mass 
star-forming regions (hereafter HSRs), such as M16 (Hill et al. 
2012; Tremblin et al. 2014), M17 (Felli et al. 1984; Chen et al. 
2021; Yin et al. 2022), W31 (Beuther et al. 2011; Gama et al. 
2016; Maity et al. 2022), W33 (Messineo et al. 2015; Tursun 
et al. 2022), and W39 (Kerton et al. 2013). At the same time, 
the infrared dust clumps of Herschel (Elia et al. 2017) also 
cover the HSR, providing us with an opportunity to validate the 
Facet-SS-3D-Clump using data from different wavelength 
bands. The dense molecular cores and clumps are defined as 
compact (~0.1 and 1 pc, respectively) and dense (>10*-10° 
H5 cm ?) structures (e.g., Williams et al. 2000; Zhang et al. 
2009; Ohashi et al. 2016; Motte et al. 2018). Here, we refer to 
the compact regions traced by CO as “clumps.” Using Facet- 
SS-3D-Clump, we obtained a catalog of clumps in the HSR. 

The structure and organization of the article are as follows: 
Section 2 introduces the observations of MWISP. Section 3 
describes the details of Facet-SS-3D-Clump that contains the 
molecular clump extraction algorithm, parameter estimation 
algorithm, completeness experiments, and SS-3D-Clump to 
verify clump candidates. In Section 4, we present a catalog 
containing 18,757 8Co clumps extracted from HSR and 
provide a brief statistical analysis of clumps in terms of their 
spatial distribution and the matching results with Herschel dust 
clumps. The conclusion is provided in Section 5. 


2. Data Introduction 


The observations of MWISP use a nine-beam superconduct- 
ing spectroscopic array receiver, working in the sideband 
separation mode and employing a fast Fourier transform 
spectrometer (Shan et al. 2012; Su et al. 2019). The lines of 
three CO isotopologues, including ?CO, ?CO, and C'%0, can 
be observed simultaneously. All observations are taken in 
position-switch On-The-Fly (OTF; see Sun et al. 2018) mode. 
The observing strategy and data reduction adopted for the 
MWISP Survey are described in detail in Su et al. (2019), and 
the details on the telescope can be found at http://www. 
radioast.nsdc.cn/mwisp.php. The CO line was observed at 
the upper sideband with a main beamwidth of ~48”, and the 
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Figure 1. The noise level of the observational data. (a) Histograms of the rms noise levels (o,;) for "CO (blue), ^cO (green) and c'80 (red). (b) Cumulative 


distributions of C,ms for three spectra. 


'3CO and C'5O lines were observed at the lower sideband with 
a main beamwidth of ~50”. Both bandwidths are 1000 MHz 
wide with 16,384 channels, resulting in velocity separations of 
about 0.159 km s^! for "CO and 0.166 kms”! for "CO and 
C50. The typical noise temperatures including the atmosphere 
at 110 GHz and 115 GHz are 140 K and 250 K, respectively. A 
detailed description of the instrument was given by Shan et al. 
(2012). 

The sky coverage of the MWISP project is divided into 
10,941 cells for the region. Each cell, being 30' x 30’, is 
scanned along Galactic longitude (/) and Galactic latitude (b) at 
least twice to reduce the fluctuation of noise. The reasons for 
choosing a cell size of 30’ x 30' are detailed in Su et al. (2019). 
The sample spacing between adjacent scans is 15". Therefore, 
the MWISP survey fully sampled the mapped area. The OTF 
raw data were regridded into 30" x 30"pixels. The data 
reduction was carried out using the GILDAS software? 
described in detail in Pety (2005). Finally, the three-dimen- 
sional (3D) FITS data cubes of each cell were made with a grid 
spacing of 30" for '7CO, ?CO, and C'^O (J = 1 — 0) lines. The 
rms distributions are presented in Figure 1. The typical rms 
noise levels of the spectra are 0.47, 0.22 and 0.21K for 
PcoJg-1-0, "CcoQ-1-0) and C¥OV=1-0), 
respectively. 

The CO integral / — b map of HSR is shown in Figure 2, 
which is integrated from visa = —5to 80km s |. The inte- 
grated intensity threshold is 2 times c,,4,, where o,,,, is the rms 
noise per velocity channel. There are several renowned star- 
forming regions in HSR. Surveying the molecular clumps in 
HSR can provide a substantial sample for studying different 
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evolutionary stages of molecular clumps. It is also a crucial step 
in understanding the early stages of stellar evolution. 


3. Generation of Clump Catalogs from MWISP 


Star-forming regions are messy and chaotic environments, 
with structures on many scales (Wurster & Rowan 2023). The 
entire region is typically referred to as a cloud, dense regions 
embedded within the cloud are clumps, while the very dense 
regions in the clumps are cores (Bergin & Tafalla 2007). There 
is general agreement of the three terms, although there is 
ambiguity among the specific definitions and divisions between 
the levels. For example, Rathborne et al. (2009) identified MCs 
and clumps using the ClumpFind algorithm, utilizing 
PcO (/—1-— 0) emission line data observed by the GRS 
(Jackson et al. 2006). Takekoshi et al. (2019) analyzed the 
statistical properties of c5og-1-0) clumps in the Cygnus 
X cluster-forming region using data from the Nobeyama 45 m 
radio telescope, with the clumps identified through the 
dendrogram algorithm (Rosolowsky et al. 2008). Liu et al. 
(2022) refer to centrally concentrated structures as clumps and 
treat dendrogram-defined leaves accordingly. 

The main characteristics of the molecular clumps are local 
intensity enhancement and different shapes. Meanwhile, they 
are embedded in lower average bulk density molecular gas 
(Blitz & Stark 1986; Lada 1992). The clumps (or cores) are 
empirically defined as the regions with concentrated, enhanced 
intensity in a data cube. It is worth noting that CO can only 
trace structures of specific density, which falls within the range 
of 107—10* cm ?. This limitation results in ^CO only revealing 
local structures within MCs. Therefore, in this study, we define 
the structures traced by CO data as clumps. 

The workflow of Facet-SS-3D-Clump for obtaining the 
catalog of molecular clumps is shown in Figure 3, in which we 
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Figure 2. The intensity maps of '*CO emission for the HSR in velocity ranges from —5 to 80 km s~". The integrated intensity threshold is 2 x gyms, where Omg is the 


rms noise per velocity channel. 


The Catalogue of Clumps 


Figure 3. Workflow of obtaining the catalog of clumps. 


utilized the dendrogram (Rosolowsky et al. 2008) and 
FacetClump (Jiang et al. 2023) algorithms to detect '^CO data. 
Subsequently, we performed cross-matching to obtain clumps 
with high confidence. Simultaneously, these high-confidence 
clumps were used as prior knowledge to train an SS-3D-Clump 
model for verifying the remaining clumps, resulting in 
confidence levels for molecular clumps. At the same time, we 


used an MGM described in Luo et al. (2022) to simultaneously 
fit overlapping molecular clumps and derive the geometric 
morphological parameters of those clumps. Finally, we 
obtained the catalog of molecular clumps, which includes the 
positions, morphological parameters, and confidence levels 
associated with the clumps. 

To avoid the difficulty of detecting clumps from the entire 
MWISP data, which is computationally expensive, we perform 
the clump extraction with a strategy that processes sub-cubes 
and global stitching (see Figure 4). A sub-cube spans 
120 pixels in the longitude and latitude directions, i.e., ~1°, 
given the MWISP data pixel size of 30". The candidates are 
rejected from the first catalog since their centers are outside the 
local region of the sub-cube (see Figure 4). Then, the clumps 
belonging to the entire region are obtained by retaining the 
clumps whose centers are in the local regions of all sub-cubes. 
As shown in Figure 4, the clumps labeled with designations a, 
b, c are in sub-cube 1, while the center of the clump c marked 
by a red dot is out of the local region of sub-cube 1. Thus, 
clump c colored with gray is removed from sub-cube 1 and 
clumps a and b colored with blue are retained in the catalog of 
sub-cube 1. In sub-cube 2, the center of clump b is out of the 
local region, while that of clump c is in the local region. The 
clump c is retained in the catalog of sub-cube 2. 


3.1. Clump Extraction Algorithm 


Dendrogram (Rosolowsky et al. 2008) implemented in 
astrodendro® is an abstraction of the changing topology of 
the isosurfaces as a function of contour level, which uses a tree 


© https: //dendrograms.readthedocs.io/en/stable / 


Research in Astronomy and Astrophysics, 24:055018 (17pp), 2024 May 


Figure 4. Scheme of processing sub-cubes and global stitching for the entire 
MWISP data. Sub-cube edges are indicated with black vertical lines, which 
span all approximately 1? in longitude (latitude is the same). The light blue 
region surrounded by vertical dotted lines marks the local region of the sub- 
cube. Clumps a and b in sub-cube 1 and clump c in sub-cube 2 are retained, 
according to whether the center of the clump falls into the local region. 


diagram to describe hierarchical structures over a range of 
scales in a two-dimensional (2D) or 3D datacube (Zhang et al. 
2021). There are two types of structures returned in the results: 
leaves, which have no substructure, and branches, which can 
split into multiple branches or leaves. There are two main 
parameters in the algorithm: Tmin and AT. Tmin is the minimum 
value to be considered in the data set. In the fiducial case, we 
adopt Tii, = 30rms, where Gyms is the local noise level. AT 
describes how significant a leaf has to be in order to be 
considered an independent entity. We adopt a fiducial value of 
AT =20,ms, which means a clump must have a peak flux 
reaching 50;ms above the noise. 

FacetClumps’ was developed by Jiang et al. (2023), which 
uses the Gaussian facet model to fit local surfaces and employs 
the extremum determination theorem of multivariate functions 
to determine clump centers. Based on the identified clump 
centers, FacetClumps clusters the regions near the centers by 
considering connectivity and minimum distance, thereby 
obtaining and detecting molecular clumps. The FacetClumps 
algorithm uses default parameters in this paper (see Appendix 
C in Jiang et al. 2023). By cross-matching the detection results 
of the two algorithms, “extremely” high-confidence clumps can 
obtained. As shown in Figure 5, drawn using the Cube 
Analysis and Rendering Tool for Astronomy (CARTA; Comrie 
et al. 2021), the green circles and blue plus signs represent the 
centroids of clumps detected by FacetClumps and dendrogram, 
respectively. When the deviations of the peak positions of two 
clumps detected by two algorithms are less than 2 pixels in all 
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Figure 5. Dendrogram and FacetClumps were applied to the same CO data; 
the green circles mark the centroids of clumps detected by FacetClumps, while 
the orange plus symbols represent the results by dendrogram. 


three directions, corresponding to deviations less than 1’ in 
Galactic longitude and latitude and less than 0.32 km s lin the 
velocity axis in the actual data, they are considered to be the 
same clump. Among them, FacetClumps and dendrogram 
detected 24,367 and 23,019 clumps in HSR, respectively, with 
10,754 clumps matched between the two. 


3.2. Verified by SS-3D-Clump 


To minimize false positives, both Rigby et al. (2019) and 
Rojas et al. (2022) employed visual inspection by independent 
reviewers, assigning reliability flags for CHIMPS clumps and 
classifying lens candidates into two catalogs using high- 
resolution imaging and spectroscopy, respectively. Addition- 
ally, to minimize the introduction of subjectivity in manual 
verification, the researcher makes every effort to control this 
uncertainty during the verification process. For example, Rojas 
et al. (2022) re-displayed the previously classified ring galaxies 
during the verification process as a consistency check (users 
should re-classify them as rings, or at least not classify them as 
lenses). Therefore, manually verifying candidates requires 
significant effort and time, especially for large amounts of data. 

Therefore, an automated verification method as a substitute 
for manual verification becomes increasingly necessary in the 
era of extensive data. Inspired by the success of supervised 
deep learning in galaxy classification (Zhu et al. 2019; Cheng 
et al. 2020; Lukic et al. 2019; He et al. 2021; Gupta et al. 
2022), we integrate limited labeled data with a deep learning 
approach. Utilizing semi-supervised deep learning, namely SS- 
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Figure 6. The structure of SS-3D-Clump is illustrated in the figure. The proposed method refers to iteratively clustering deep features and using cluster assignments as 
pseudo-labels to learn the parameters of the 3D Convolutional Neural Network (CNN) and classifier networks. 


3D-Clump (Luo et al. 2024), we verify the clump candidates. 
As shown in Figure 6, the deep features of the clump candidate 
are extracted by the feature extraction part of SS-3D-Clump 
(Luo et al. 2024). Utilizing the deep features as a foundation, 
the Constrained-KMeans (Basu et al. 2002) algorithm effi- 
ciently generates pseudo-labels for candidates using small 
labeled samples as seeds. The SS-3D-Clump classifier also 
assigns predict-labels to the candidates based on their deep 
features. The difference between both labels could be used to 
optimize the parameters of the SS-3D-Clump model. SS-3D- 
Clump iteratively groups the deep features with a standard 
clustering algorithm and uses the subsequent assignments as 
supervision to update the network weights of SS-3D-Clump. 
Figure 7 illustrates the curve of the Normalized Mutual 
Information (NMI) changing with the epochs during the 
training process of SS-3D-Clump. NMI is a metric used to 
assess the similarity of two data sets or the performance of 
clustering algorithms. The range of NMI typically falls between 
O and 1, where 1 indicates complete similarity between two 
clustering results, and O signifies no resemblance. During the 
training process, NMI curves can be utilized to observe how the 
algorithm's performance changes with increasing training 
epochs, facilitating the optimization of the model training 
process. Measuring NMI between the clusters at epoch t— 1 
and f provides insights into the actual stability of SS-3D- 
Clump. From Figure 7, it can be observed that, by the 15th 
epoch of training, the NMI increases to around 0.97 and 
remains stable. This indicates that SS-3D-Clump is experien- 
cing fewer reassignments, and the clusters are stabilizing. 
Since SS-3D-Clump belongs to semi-supervised deep 
learning, the verification of molecular clumps is completed 
simultaneously with the model training. Utilizing the stabilized 
SS-3D-Clump after training, we conducted verification on the 
candidates detected by FacetClumps, excluding those with 
confidence levels below 0.8. In the end, SS-3D-Clump 
excluded 5792 out of the 24,367 molecular clump candidates 
detected by FacetClumps, resulting in a catalog containing 
18,575 molecular clumps. The histograms of peak intensities of 
the molecular clumps are illustrated in Figure 8. The green 
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Figure 7. Evolution of cluster reassignments (NMI) in each clustering iteration 
of SS-3D-Clump. 


histogram depicts the peak value distribution of clumps 
retained by SS-3D-Clump, while the red one shows the cases 
excluded by SS-3D-Clump. It can be inferred from Figure 8 
that the distribution of retained clumps has a sharp rise spike 
with position of about 1.5 K. After passing the peak position, 
the distribution drops rapidly and slowly descends. 
Appendix A presents a subset of the results obtained through 
SS-3D-Clump verification, which contains integrated maps of 
molecular clumps on the / — b plane along with their associated 
confidence levels. 


3.3. Morphological Parameter Estimation 


For blended source pairs, the Gaussian fitting results in a 
larger size estimate, whereas in cases where two or multiple 
sources are resolved by the detection algorithm, simultaneous 
fitting of multiple components enables better separation 
(Molinari et al. 2016). Therefore, we use the MGM? described 
in Luo et al. (2022) to simultaneously fit overlapping molecular 
clumps and derive the geometric morphological parameters of 
those clumps. The integrated map shown in Figure 9 represents 
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Figure 8. Histogram illustrating the distribution of peak values for molecular 
clump candidates detected by FacetClumps and verified by SS-3D-Clump. The 
green bin markings represent the distribution of retained clumps, while the red 
ones represent those excluded by SS-3D-Clump. 


a local region in HSR with an integrated velocity range of 
25.6-28.] kms !. For the molecular clumps detected by 
FacetClumps, morphological parameters are estimated using 
MGM. The green dots signify the centroids of clumps, and the 
white ellipses are drawn based on the clumps’ major axis, 
minor axis, and rotation angle. 


3.4. Completeness 


The limitation of the telescope sensitivity causes a low 
quality clump to be ignored. Other indicators of the algorithm 
are the completeness and the detection rate above the 
limitation. The “completeness limit” here refers to the total 
flux or mass above which a clump can be detected at a certain 
level with an algorithm. The smaller and weaker molecular 
clumps are, the less likely they are to be detected. To quantify 
the degree of completeness of the extracted clumps, we carried 
out an extensive set of synthesis data experiments by injecting 
simulated clumps into the “CO(J=1-—0) data in HSR. 
Completeness experiments for '*CO were conducted by 
partitioning the entire HSR into 16 subregions, collectively 
covering the entire HSR. 

Via multiple injections, we injected 36,000 simulated clumps 
modeled as 3D ellipsoidal Gaussians for each subregion, with 
their axis size and intensities randomly distributed within 
specific ranges. The peak intensity value of those simulated 
clumps takes values from 0.7 to 15, while the size of the clump 
in the velocity axis takes values from 2 to 5 and the size in the 
Galactic longitude and latitude axes takes values from 1.5 to 4. 
In this way, we were able to assess the capability of recovering 
a statistically comparable population of clumps from the 
designated region. These simulated clumps were randomly 
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Figure 9. The figure displays an integrated map of a local region in HSR with 
an integrated velocity range of 25.6-28.1 kms’. Morphological parameters 
for the molecular clumps detected by FacetClumps are estimated using MGM. 
The green dots indicate the centroids of clumps, and the white ellipses are 
drawn based on the major axis, minor axis, and rotation angle of the clumps. 


spread on the map, with the only constraint being to avoid 
overlap in their positions. 

The synthesis data were processed with FacetClumps, 
adopting the same parameters as used for real data of HSR, 
and the outputs were compared with the truth table of the 
simulated clumps. The number of clumps within each total flux 
interval was counted, and the completeness was obtained. An 
example of the recovery fraction as a function of the integrated 
flux density of the simulated clumps in different latitude 
intervals in HSR is shown in Figure 10. From Figure 10, it can 
be observed that the completeness curves decrease significantly 
when the flux of clumps is from 40 to 20 Kkms ! within 
different latitude intervals. When the completeness reaches 
90% for the entire region, the total flux of clumps is 
approximately 37 K kms '. In Figure 11, the dashed line 
represents the estimated completeness limit as a function of 
Galactic latitude. At the same time, the histogram shows the 
distribution of missed simulated clumps in Galactic latitude 
during the synthetic data experiment. Around 0^, specifically in 
—1?5« b < 075, the dashed line exhibits a noticeable increase, 
corresponding to a peak in the histogram at this position. It is 
explained by the brighter emission at lower latitudes, making 
detecting fainter objects a more challenging task. 

The completeness limits reported in Figure 11 should be seen 
as conservative because they are determined by spreading the 
synthetic clumps randomly over each entire region. However, 
the background has non-uniform characteristics in each region, 
decreasing toward the north and south Galactic directions as the 
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extracted catalog with statistically the same sizes. 
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Figure 11. Ninety percent completeness limits in flux density for a population 
of clumps with the same distribution of sizes as the one extracted by 
FacetClumps as a function of Galactic latitude. The significant increase in the 
completeness limit near low Galactic latitudes is due to the brighter background 
emission in HSR. 


Galactic latitude increases. In Figure 12, a typical example is 
shown, with the upper panel showing the integrated CO map 
of MWISP centered at / — 17°. Superimposed are the extracted 
clumps with integrated fluxes above (yellow crosses) and below 
(blue crosses). The lower panel of Figure 12 further depicts the 
latitude distribution of the two groups of clumps, also signified 
with solid and dashed lines for sources above and below the 
confusion limit respectively. The two groups of clumps exhibit 
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Figure 12. The spatial distribution of molecular clumps above and below the 
completeness limit and their distribution in Galactic latitude. Upper panel: the 
integrated "CO map of MWISP at |= 17°. Superimposed are the clumps 
detected by FacetClumps. The yellow crosses indicate the clumps with fluxes 
above the completeness limit, while the blue crosses indicate the clumps below 
the completeness limit. Lower panel: Histograms of latitude distributions for 
clumps, above (solid line) and below (dashed line) the completeness limit. 


subtle differences in their spatial distribution, with clumps 
brighter than the completeness limit mainly concentrated in the 
Galactic latitude range of — 1? to 076. In contrast, the distribution 
of the other group of clumps is more uniform. 
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Figure 13. The one-to-one matching results of Herschel 350 jum dust clump and CO clump. The upper-left subplot displays the channel map with the maximum 
intensity of the '^CO clump, with the corresponding Herschel dust clump located in the upper-right subplot. The lower subplot illustrates the average spectrum of the 
13CO clump, and the averaging range is delineated by the black circle in the upper-left subplot. 


Column Name 


ID 

GLON 
GLAT 
Velocity 
FWHMA 
FWHMB 
FWHMV 
PA 

Peak 

Flux 

S/N 

Fitting Error 
Confidence 


Table 1 
Description of Columns in the Catalog 
Units Symbol Explanations 
Designation of the clump (MWISPLLL.IIIJ-B.bbb--V V V.v) 
deg L Galactic longitude 
deg b Galactic latitude 
kms! VLSR the central radial velocity with respect to LSR 
arcsec Su FWHM of the clump along axis a of the elliptical Gaussian as determined by fitting engine 
arcsec Sin FWHM of the clump along axis a of the elliptical Gaussian 
kms! S, FWHM of the clump in the velocity direction 
deg 0 Position angle 
K Tpeak The Peak value of clump 
Kkms^! The total integrated flux of clump 


wW 
r Signal-to-noise ratio obtained by dividing Peak by the residual rms 

E K E is calculated by the root mean square error between the fitting and the actual value 
a The confidence associated with the clumps verified by SS-3D-Clump 
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Figure 14. The figure illustrates results similar to Figure 13, with the distinction that it showcases a Herschel dust clump matched to multiple '*CO clumps. 


4. Results 
4.1. Detailed Contents of Clump Catalog 


A total of 18,575 "CO clumps are extracted in HSR. 
Table B1 in Appendix B presents a part of the clump catalog, 
while the full catalog can be obtained online (https:/ /www. 
scidb.cn/en/s/qEfe2m). In the table, columns (2) and (3) are 
the positions of the clumps in galactic coordinates, and column 
(4) is the central radial velocity with respect to the local standard 
of rest (LSR). Columns (5) and (6) represent the major and minor 
sizes of the fitted ellipse on the sky plane, respectively. Column 
(7) represents the FWHM size in the velocity direction. Column 
(8) 1s the tilt angle on the sky plane, calculated from the north 
Galactic pole to the major axis of the clump in the counter- 
clockwise direction. Columns (9) and (10) list the clump's peak 
intensity and total integrated intensity value, respectively. 
Column (11) represents the peak signal to noise ratio (S/N) of 
the clump. Column (12) is the fitting error, which is calculated by 
the rms error between the fitting value and the actual value. 
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Column (13) is the confidence levels associated with the clumps, 
which are verified by SS-3D-Clump. Names, units, symbols, and 
detailed descriptions of the columns in the molecular clump 
catalog are provided in Table 1. 


4.2. Cross Match with Hi-GAL Sources 


Herschel infrared Galactic Plane Survey (Hi-GAL) is a large- 
scale survey of the Galactic plane (Molinari et al. 2010), 
performed with Herschel in five infrared continuum bands 
between 70 and 500 jum using the Herschel Space Observatory 
(Pilbratt et al. 2010). The first public release of the high-quality 
products from the Hi-GAL survey was presented by Molinari 
et al. (2016), in which ~10° compact sources are identified in 
each band. Elia et al. (2017) merged the sources of the five 
bands mentioned above to obtain a band-merged catalog 
containing 100,922 sources with a regular SED, 24,584 of 
which show a 70,un counterpart and are considered 
protostellar. 


Research in Astronomy and Astrophysics, 24:055018 (17pp), 2024 May 


T. 1e3 


L—1 BPCO Clump 
C] Hi-GAL Clump 


1.5 


1:2 


0.9 


0.6 


Number of Clumps 


0.3 


0.0 Poe 


3 -2 -1 0 1 2 3 


Galactic Latitude [deg] 


Figure 15. Distribution of molecular clumps along Galactic latitude in HSR. 
The blue histogram represents the distribution of molecular clumps detected by 
Facet-SS-3D-Clump, while the orange histogram represents the distribution of 
sources in the Herschel survey. 


Within the Galactic longitude range 10? « / « 20°, there are 
8517 band-merged sources, and through matching with CO 
clumps, 7154 sources have been identified as corresponding to 
'3CO clumps. Among them, 2200 sources are one-to-one 
matches with ^CO clumps, while 4954 sources match multiple 
?cO clumps at different velocities. As shown in Figure 13, the 
figure illustrates the one-to-one matching results of the 
Herschel 350 jum dust clump and '*CO clump. The upper-left 
subplot of Figure 13 shows the channel map with the maximum 
intensity of the SCO clump and the corresponding Herschel 
dust clump is in the upper-right subplot. The lower subplot of 
Figure 13 represents the average spectrum of the SCO clump, 
with the averaging range indicated by the black circle in the 
upper-left subplot. Looking at the blue contours of CO in 
Figure 13 and the overlay of these contours onto the dust 
clumps, the positions of the two correspond very well. As 
shown in Figure 14, similar to Figure 13, the only difference is 
that the displayed results depict a Herschel dust clump matched 
to multiple "CO clumps. The figure reveals a similarity 
between the CO data and Herschel dust emission. In the one- 
to-multi matching results, the average spectrum of CO data at 
the positions of Herschel dust clumps exhibits multiple peaks, 
which can explain why the number of CO clumps is greater 
than that of Herschel dust clumps. The upper-left subplot of 
Figure 14 displays the channel map of the maximum intensity 
in the average spectrum. 


4.3. Spatial Distribution 


A comparison with Hi-GAL dust clumps (Elia et al. 2017), 
covering the inner Milky Way in the longitude range 
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Figure 16. The distribution of CO clumps along the Galactic longitude, 
where the dashed line signifies the star-forming region in the HSR. The figure 
reveals a local clustering of CO clumps near the longitude of the star-forming 
region, indicating the consistency of clumps detected by Facet-SS-3D-Clump. 


10*?«1«20?, shows substantial similarities in the source 
count distributions. The latitude distribution of the HSR clumps 
and the Hi-GAL compact sources is reported in Figure 15. The 
blue histogram represents the distribution of molecular clumps 
obtained by Facet-SS-3D-Clump, while the orange represents 
the distribution of dust clumps in Hi-GAL (Elia et al. 2017). 
Both histograms peak at slightly negative values, and the 
median value for HSR clump latitude is ~ —0°118 below the 
nominal midplane. Also, the mean Galactic latitude for clumps 
within |b| < 2? is significantly below the midplane: 
b — —07110. From the histograms, it can also be observed 
that the distribution of molecular clumps exhibits an asymme- 
trical pattern. However, we cannot rule out the asymmetric 
distribution with respect to the central plane in the studied field. 
As Molinari et al. (2016) pointed out, this may be attributed to 
an incorrect assumption about the vertical position of the Sun 
in the Milky Way. 

The longitude distribution of SCO clumps is shown in 
Figure 16, with the dashed line representing the star-forming 
region in HSR. From the figure, it can be observed that there 
is a local clustering of 13CO clumps near the longitude of the 
star-forming region, which also reflects the consistency of 
clumps obtained by Facet-SS-3D-Clump. The top panel of 
Figure 17 shows the positions of clumps in the / — b plane. 
The bottom of Figure 17 presents the centers of clumps in the 
l—v plane. The yellow, red, green, and purple lines 
correspond to the spiral arms of Aql Rift, Sagittarius, 
Scutum, and Norma, respectively. Dashed lines indicate 
the nearness of the arms, while solid lines represent their 
farness. 
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Figure 17. The distribution of clumps in both spatial and velocity dimensions. Top: the positions of clumps in the / — b plane. The green dots indicate centroids of 
clumps. Bottom: the positions of the clumps are shown on the / — v plane. The yellow, red, green, and purple lines represent the Aql Rift, Sagittarius, Scutum, and 
Norm's arms, respectively. Dashed lines indicate the nearness of the arms, while solid lines represent their farness. 
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5. Conclusion 


We employed Facet-SS-3D-Clump to identify and verify 
molecular clumps in the CO data of MWISP, covering the 
longitude range 10°</<20° and a latitude strip of 
|b| < 5? 25. Through the analysis of the results, we conclude 
that: (1) We obtained a catalog containing 18,575 '?CO clumps 
from the MWISP data using Facet-SS-3D-Clump. (2) The flux 
completeness limits showed that the catalog is 90% complete 
above 37 Kkms ’. (3) There are 82.3% of the dust clumps 
corresponding to '*CO clumps from matching results with the 
Herschel infrared dust clumps. The Facet-SS-3D-Clump will be 
applied to survey SCO clumps across the entire MWISP data. 
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Appendix A 
Examples of Clumps verified by SS-3D-Clump 


To provide a demonstration of the efficacy of SS-3D-Clump 
in verifying molecular clump candidates, a subset of the 
samples were visualized. Figure A1 shows molecular clumps 
on the /-b plane with confidence levels obtained through SS- 
3D-Clump verification. 
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Figure A1. The integrated maps on the / — b plane of the molecular clumps are presented. In the figures, white numbers indicate the confidence assigned by SS-3D- 
Clump for the current clump, blue dots represent the centroids of the clumps, and white ellipses depict the results of the 3D Gaussian modeling performed by MGM 
for the clump. The white ellipses mark the positions of the Gaussian ellipsoidal FWHM. 
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Appendix B 


Partial Molecular Clumps in the Catalog 


Table B1 presents a partial catalog of the clumps, providing 
information about their positions, sizes, intensities, and fluxes. 
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Table B1 
Partial Catalog of Molecular Clumps 
ID l b VLSR Su Sn S, 0 Tpeak W E a 
MWSIPO12.227-00.031+032.3 12.227 —0.031 32.3 136 101 1.18 7 3 65 1.32 0.99 
MWSIPO012.091-00.069+013.1 12.091 —0.069 13.1 148 101 1 180 22, 42.7 0.69 0.99 
MWSIPO11.845-00.100+015.6 11.845 —0.1 15.62 160 118 1.63 180 6.4 261.8 1.34 0.99 
MWSIPO015.583--00.269--035.5 15.583 0.269 35.5 142 93 0.43 38 1.1 8.5 0.48 0.57 
MWSIP012.505-00.2364-004.3 12.505 —0.236 4.31 168 124 0.85 90 12 29.6 0.57 0.85 
MWSIPO012.033-00.1214-022.2 12.033 —0.121 22.16 211 99 0.89 27 0.9 21.3 0.95 0.96 
MWSIPO012.200-00.0394-050.5 12.2 —0.039 50.47 135 124 1.9 136 6.4 268.8 1.22 0.97 
MWSIPO011.879-00.0874-050.3 11.879 —0.087 50.27 121 88 1.32 125 1.4 25.6 0.59 0.94 
MWSIPO015.4674-00.2724-029.0 15.467 0.272 28.98 124 100 0.98 180 1 16.4 0.53 0.86 
MWSIPO012.010-00.1814-038.9 12.01 —0.181 38.94 262 176 1.43 132 3.8 338 1.46 0.99 
MWSIPO012.031-00.0324-110.7 12.031 —0.032 110.71 125 105 0.87 48 6.4 96.7 1.32 0.98 
MWSIPO12.031-00.250+117.4 12.031 —0.25 117.39 113 103 1.37 118 2.1 45.5 0.70 0.99 
MWSIPO11.776-00.165+044.4 11.776 —0.165 44.4 134 61 0.69 35 1.8 13.5 1.12 0.98 
MWSIP011.848--00.005--012.3 11.848 0.005 12.32 198 116 1.09 129 1.6 52:5 0.74 0.97 
MWSIP012.093-00.1024-113.1 12.093 —0.102 113.09 188 104 0.98 129 1.7 42.7 0.84 0.99 
MWSIPO017.513-01.1144-043.3 17.513 —1.114 43.3 130 66 0.33 140 1.0 3.6 0.57 0.46 
MWSIPO11.966+00.050+018.9 11.966 0.05 18.91 151 93 1.16 37 2.6 55,5 0.77 0.99 
MWSIP012.058-00.1424-019.6 12.058 —0.142 19.58 135 91 0.71 20 12 14 0.54 0.58 
MWSIPO012.060-00.1504-000.8 12.06 —0.15 0.77 135 113 0.76 180 2.1 32.4 0.79 0.97 
MWSIP012.133--00.074--017.8 12.133 0.074 17.81 134 128 1.11 90 24 60.3 0.97 0.99 
MWSIPO11.766-00.068+060.4 11.766 —0.068 60.38 134 88 0.77 180 7.1 86.1 1.68 0.98 
MWSIPO11.864-00.231+046.4 11.864 —0.231 46.39 129 93 0.78 54 1.3 15.6 0.63 0.96 
MWSIP012.190--00.036--048.5 12.19 0.036 48.49 174 128 1.26 18 3.5 131.9 1.13 1.00 
MWSIP012.047--00.006--061.5 12.047 0.006 61.47 75 59 0.72 126 3.2 13.2 0.70 0.96 
MWSIP012.129+00.059+112.4 12.129 0.059 112.38 119 109 1.27 90 1.6 36.2 0.64 0.97 
MWSIPO12.213+00.014+031.4 12.213 0.014 31.43 150 136 1.12 7 3.5 105.6 0.95 0.99 
MWSIPO012.174--00.227--029.2 12.174 0.227 29.18 109 80 0.83 11 1.3 12.9 0.76 0.97 
MWSIP013.909-04.638+011.4 13.909 —4.638 11.36 105 71 0.37 1 1.6 59 0.55 0.82 
MWSIP012.016-00.210-003.5 12.016 —0.21 -3.48 137 96 1.79 100 6.5 203 1.05 0.97 
MWSIP018.812+00.272+018.5 18.812 0.272 18.5 140 90 1.03 90 3 51.8 0.99 0.99 
MWSIP013.731-00.799+019.3 13.731 —0.799 19.28 156 111 0.91 120 1.8 37.4 0.65 0.99 
MWSIP012.201-00.117+027.1 12.201 —0.117 27.06 147 140 2.14 1 8.9 521.5 1.47 0.98 
MWSIP011.884-00.131+018.8 11.884 —0.131 18.77 115 107 1.39 124 1.8 42.1 0.80 0.98 
MWSIP011.964-00.136+032.6 11.964 —0.136 32.56 153 145 0.94 130 4 112.3 0.80 0.98 
MWSIP014.313-00.825+013.4 14.313 —0.825 13.44 137 87 0.88 32 2.1 28.7 0.79 0.99 
MWSIP018.358-00.006+125.3 18.358 —0.006 125.33 113 93 0.61 113 1.5 13 0.55 0.87 
MWSIP012.231+00.193+015.6 12.231 0.193 15.62 187 162 1.29 106 2.9 152.2 0.89 1.00 
MWSIP011.970-00.209+120.2 11.97 —0.209 120.15 110 86 0.96 33 2.1 25.3 0.64 0.98 
MWSIP012.118+00.033+095.7 12.118 0.033 95.69 188 132 1.14 151 21 79.9 0.72 0.99 
MWSIP011.835-00.191+001.6 11.835 —0.191 1.65 150 71 0.84 45 1.6 18.5 0.66 0.99 
MWSIPO11.904-00.251-003.9 11.904 —0.251 -3.91 133 98 1.39 1 2 46.9 0.68 0.99 
MWSIP016.647-01.2614-016.8 16.647 —1.261 16.76 148 116 0.86 109 3.4 66.4 0.86 0.99 
MWSIPO12.106-00.027+111.6 12.106 —0.027 111.58 138 106 1.25 124 3T 89.5 1.01 1.00 
MWSIPO011.9674-00.1194-097.8 11.967 0.119 97.8 180 98 0.81 43 2.1 39.4 0.67 0.99 
MWSIP011.789-00.013+037.4 11.789 —0.013 37.35 174 112 1.02 133 1.3 332 1.34 0.97 
MWSIPO014.400-00.8244-021.4 14.4 —0.824 21.44 141 115 1.2 180 7 181.8 1.93 0.99 
MWSIP014.277-00.9444-021.8 14.277 —0.944 21.77 132 125 0.79 148 5.7 98.8 0.95 0.97 
MWSIPO011.870-00.1104-004.1 11.87 —0.11 4.07 161 124 1.26 90 2.2 71.9 0.70 1.00 
MWSIP012.194-00.024+014.1 12.194 —0.024 14.12 153 94 1.14 114 22 49.5 0.75 1.00 
MWSIPO015.2874-00.2574-028.1 15.287 0.257 28.15 163 124 0.56 155 1.4 21.4 0.49 0.82 
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Table B1 
(Continued) 

ID l b VLSR Sm Sm S, 0 Tpeak Ww E a 

MWSIPO010.072--01.092--027.7 10.072 1.092 27.69 111 89 0.97 142 1.3 17.1 0.51 0.84 
MWSIPO019.224--00.385--020.2 19.224 0.385 20.22 152 95 0.94 125 3.6 65.6 1.45 0.98 
MWSIP012.150-00.0244-047.0 12.15 —0.024 46.97 168 126 1.09 154 3.9 120.4 1.20 1.00 
MWSIP011.920--00.068--095.8 11.92 0.068 95.77 181 128 1.11 180 2 69.9 0.66 0.99 
MWSIPO011.915--00.012--034.2 11.915 0.012 34.25 192 124 1.14 147 3 106.5 1.09 0.99 
MWSIPO012.123--00.226--049.3 12.123 0.226 49.33 124 116 1.41 180 2.2 60.5 0.64 0.99 
MWSIP016.324-00.2184-025.7 16.324 —0.218 25.67 115 63 0.62 133 1.7 10.4 0.52 0.54 
MWSIP014.495-00.8734-016.0 14.495 —0.873 16.02 204 97 0.96 123 1.8 46.7 0.71 0.98 
MWSIP012.253-00.2194-047.5 12.253 —0.219 47.49 104 97 0.72 90 2:3 22.2 0.60 0.97 
MWSIP016.573--00.042--045.3 16.573 0.042 45.27 133 67 0.65 156 1.4 10.7 0.53 0.84 
MWSIPO017.647-00.7964-057.6 17.647 —0.796 57.63 98 66 0.66 8 1.5 8.6 0.57 0.57 
MWSIPO012.231--00.105--121.0 12.231 0.105 121.05 110 68 0.77 0 1.2 9.2 0.55 0.71 
MWSIPO012.207--00.169--049.9 12.207 0.169 49.91 106 94 0.83 103 2:3 26.1 0.65 0.98 
MWSIP015.841--00.936--028.3 15.841 0.936 28.32 93 68 0.62 25 1.3 6.9 0.55 0.80 
MWSIP013.732-00.8964-021.0 13.732 —0.896 21.01 178 142 1.3 16 3.9 169.7 0.82 0.99 
MWSIPO010.779-00.1394-020.1 10.779 —0.139 20.11 130 79 0.8 13 1.1 12.2 0.54 0.82 
MWSIP012.241-00.1484-014.4 12.241 —0.148 14.36 124 111 0.82 141 7 106 0.91 0.95 
MWSIP018.287--00.129--033.6 18.287 0.129 33.62 94 76 0.73 45 1.6 10.7 0.51 0.71 
MWSIP011.798--00.076--015.8 11.798 0.076 15.82 147 108 0.94 35 4 80.4 0.98 0.99 
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